A new two-faced scalar solution and cosmological SUSY breaking. 
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We propose a possible new way to resolve the long standing problem of strong supersymmetry 
breaking coexisting with a small cosmological constant. We consider a scalar component of a 
minimally coupled N=l supermultiplet in a general Friedmann-Robertson- Walker (FRW) expanding 



universe. We argue that a tiny term, proportional to H 



10- 



in Plank's units, appearing in the 



field equations due to this expansion will provide both, the small vacuum energy and the heavy mass 
of the scalar supersymmetric partner. We present a non-perturbative solution for the scalar field 
with an unusual dual-frequency behavior. This solution has two characteristic mass scales related 
to the Hubble parameter as H 1 ^ 4 and H 1 ^ 2 measured in Plank's units. 



PACS numbers: 98.80.Cq, ll.25.-w, 04.65.+e 



I. INTRODUCTION 

A discrepancy between the vacuum energy ~ Mp pre- 
dicted by quantum field theory and the observed tiny 
cosmological constant A ~ 10~~ 122 Mp is one of the ma- 
jor problems in reconciling quantum field theory with 
cosmology. Supersymmetry could partially resolve this 
problem because the contributions to the vacuum energy 
from fermionic and bosonic loops cancel each other giv- 
ing zero vacuum energy and a zero cosmological constant. 
The unsatisfactory feature of this scenario is the fact that 
supersymmetry should be broken at a rather large scale 
(compared to the cosmological constant) to give the ob- 
served shifts in the masses of superpartners. 

We propose a possible new way to resolve this long 
standing problem of strong supersymmetry breaking co- 
existing with a small cosmological constant. We consider 
a scalar component of minimally coupled N=l supermul- 
tiplet in a general Friedmann-Robertson- Walker (FRW) 
expanding universe. We argue that a tiny term, propor- 
tional to H 2 ~ 10~ 122 in Plank's units, appearing in the 
field equations due to this expansion will influence both 
the vacuum energy and the masses of supersymmetric 
partners in a very surprising way. We were able to find 
an unusual non-perturbative solution for the scalar field 
with the double periodic behavior of non-linear oscilla- 
tor. This "two-faced" behavior provides a singular scalar 
solution with two very different frequencies correspond- 
ing to mass scales of ~ 2 • 10 _3 eV and ~ 300GeV. It is 
important to emphasize that these two mass scales are 
not related to regular spontaneous symmetry breaking 
mechanism. 



A straightforward estimation shows that the size of 
the universe changes during the time of the virtual loop 
fluctuations, ti oop ~ tp where tp is the Plank time (tp ~ 
lCP 43 s). The Plank scale is the natural choice because 
it represents the most energetic part of the fluctuation 
spectra. During this time the relative change of the size 
of the universe Rq will be proportional to ~ 1CT 



The change in the scale factor a during this time can 
be estimated using the value of the Hubble parameter 
H = | . For observed H ~ 70 km/s/Mpc - 2.3 • lO -1 ^ -1 
the change in the scale factor A (a 2 ) = (atp) 2 will be 
proportional to: 



A(a 2 ) 



H 2 {t P ) 2 -6-10" 



This means that the universe will change over the char- 
acteristic virtual loop time (Plank time tp) and that the 
interactions of virtual particles in vacuum loops should 
reflect this difference This change in the virtual parti- 
cles at the beginning and at the end of a vacuum loop 
will lead to the effective "decompensation" on a scale of 
~ 1CP 122 . This amazing coincidence between tiny value 



of 



and the ratio of ^i" ac may be related to the 



statement made by P. Dirac about the " non-random" na- 
ture of the coincidence of the large dimensionless values 
in physics 

The regular SUSY models in a flat " static" space have 
a zero vacuum energy due to the cancelation of bosonic 
and fermionic contributions However, the expansion 
of the universe could break this balance and create a 
non-zero vacuum energy density. In the expanding uni- 
verse the fermionic and bosonic vacuum fluctuations be- 
have differently and the effects on the energy density will 
be of the order of p ~ 10~ 122 Mp. The estimations of 
the scalar component contribution to the energy density 
fluctuation in the expanding universe can be found in 
multiple papers: Recently, similar result for super- 

symmetric model was derived in Q. As it follows from 
these papers the scalar energy density is p s ~ H 2 with 
Mj, « 10~ 30 Af P « 3-lCT 3 eU. The fermionic contribution 
was considered in Q ■ Following this paper the fermionic 
contribution is proportional to pf ~ Rq 3 where Ro is a 
horizon size. For the current epoch this estimation gives 
a very small value pf ps (10 _22 ) 4 eV 4 . 

In this paper we will investigate the possibility of a 
novel approach to supersymmetry breaking that will pro- 
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vide relatively high fermion-scalar mass splitting (10 2 — 
10 3 GeV) while keeping the value of the vacuum en- 
ergy density less than (10~ 2 ) 4 eV 4 . We propose an al- 
ternative approach based on a non-perturbative solution 
for scalar field in the expanding Friedmann-Robertson- 
Walker background. 



II. NON-PERTURBATIVE SCALAR SOLUTION 

General N=l supersymmetric Lagrangian of chiral su- 
permultiplets (see [2J and references therein) has a form: 



L = d^jd^i - \m,ij4>j + \ijk<f>j<f>k\ u 
i - ,,„ , i 
2< 



(1) 



where is a set of scalar components of chiral multiplets 
and ipi are Majorana spinors from the same multiplets. 
We will concentrate on a single massless scalar field with 
a quartic self-interaction in a simplest form: 

L = d»(f>d»<f>-\<l> i . 

For the expanding universe with a Hubble parameter 
H = — the equation of motion for the field <f> has a form: 

• dV 
^ + 3Hcj ) + — = 0. 

0(j) 

With the variable transformation: 

u(t) 



<P(t) 



a(t) 3 / 2 ' 

the equation of motion becomes: 





fa 




fa 


-1 




\ 





u + a 3 / 2 V'(u/a 3 / 2 ) = 0, (2) 



where an overdot indicates a derivative with respect to t 
and a prime is a derivative with respect to the field <j>. The 
second, "effective mass" linear term: 



3 

2 V2 J 



be estimated using Friedmann-Robertson- Walker equa- 
tions [1, with the metric ds 2 = —dt 2 + a 2 (t)dx 2 



4-7T 



G{ P + 3p): H 2 



-Gp. 



The linear "mass" term related to the deceleration pa- 
rameter q = sr. The current value of the decelera- 

^ aH A 

tion parameter is qo ~ —0.5 and the mass term can be 
estimated as -§ (\H 2 + f ) ps -1.5 • H 2 ps -H 2 . In the 
early universe the value for the deceleration parameter 
was probably quite different. In the case of hot relativis- 
tic gas dominance, p = p/3, and the mass term will flip 
sign as — | (^H 2 + -) w jH 2 . We will concentrate on 



the present moment, in this case the approximate equa- 
tion will take form: 



il - H 2 u + 4ru 3 
a 6 



0. 



(3) 



for simplicity we set the numerical factor ~ 0(1) in front 
of the second term to 1. This non-linear oscillatory equa- 
tion has two non-perturbative solutions: 



Ml 2 = ±i" 



• sn 



n(t — c 2 ); —1 — 



IP 

K 2 



(4) 



where sn 



nt: —1 



fi- 



ll], parameter n is 



is the Jacobi elliptic function 



Aci 



i 3 H 2 + y/a^H 4 + 2a i Xc 1 ' 



(5) 



and ci and ci are free parameters. The parameter Ci 
is an arbitrary time-shift that we will set to zero in all 
further expressions. If the mass term changes sign, the 
equation 



u - 



H 2 u 



A 



a 



3«* = 0. 



will have solutions of the form 
' -ci 



"1,2 = T«" 



• sn 



nt: —1 



H 2 



Act 



where k - y - a s H 2 +Va H H i +2a s Xci ' 

We set the value of the self-interaction constant to A s» 
10 -12 , the choice of the stronger self-interaction constant 
A fa 10 -2 is also possible, see @. This specific value of 
A will affect only the value of the free parameter c\ and 
is not crucial for the general results of this paper. The 
current value of the scale factor a can be set to a = 1, 
other possible choices discussed in The value of 

parameter c\ will be considered later. 

The Jacobi elliptic function sn[r, k] is one of the 12 Ja- 
cobi functions of this kind, [Tlj . The second argument, 
k is the elliptical modulus related to parameter m as 
k 2 = m. In our solution m = — 1 — S— with k defined in 
eq. (J5J) . The presence of a tiny second term in m pro- 
portional to H 2 will lead to very important consequences 
discussed later. The main feature of the Jacobi elliptic 
functions is the existence of two distinct complex non- 
collinear oscillatory periods. The two periods of elliptic 
sn[r, m] are defined as follows 

sn[r + ApK + ilnK' , to] = sn[r, to], 

where p and n are integer numbers, and K(m) is the 
complete elliptic integral of the first kind with K'(m) = 
K (1 — to). The real physical time t is related to a di- 
mensionless complex parameter r as r = nt. It means 
that the direction of r evolution with physical time t on 
a complex plane will be defined by k. 
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The non-linear equation of the type (|3]) for particle 
physics was previously studied by many authors, includ- 
ing M. Fraska in [T3| - In these papers, however, the sec- 
ond term of equation ([3]) proportional to H 2 was missing 
and the value of parameter m = — 1 . 



III. DOUBLE-PERIODIC SOLUTIONS 

The double-periodic nature of solution ((4]) with H 2 
dependence opens a very interesting possibility of exis- 
tence of two distinct non-linear oscillatory processes in a 
single two-faced solution for the scalar field <f> with drasti- 
cally different effective frequencies. We will show that the 
heavy component can have an effective mass of the order 
~ 300 — 400 GeV and the light component's mass can be 
of the order of ~ (2 — 3) ■ 10 _3 eV for the current epoch. 
These numbers will depend on the cosmological time and 
the corresponding value of the Hubble parameter. Inter- 
estingly, this dependence is very soft and smooth for the 
time period between the end of inflation and current time. 
During the inflation and reheating both scales will follow 
and, possibly, influence the metric evolution. 
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FIG. 2: The heavy M\ mass component of the scalar field 
eq.© for the Hubble value H = 10" 6 4 P . The dotted line 

corresponds to the set of c r and a values for 
Mi w 3 • 10 2 GeV and the dashed line indicates c r and d 
region for the light mass M2 ~ 2.8 • 10 _3 eV. 
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FIG. 1: The periods T{ and define the physical periods 
along kI direction 

The periods T\ = 4J ^(") and T 2 = l2K ( 1 ~ mS > are com _ 

plex valued and have different directions on a complex 
plane. The real physical periods T\i and T 2 ! and the cor- 
responding frequencies of the non-linear oscillator Q are 
not directly defined by values of Ti and T 2 but are re- 
lated to them through the parallelogram projections on 
a complex plane (see Figure 1). With the evolution of 
the real physical time the dimensionless argument r of 
the Jacobi sn[r, m] function evolves in the direction of k. 

The periodic structure of Jacobi sn function is plotted 
in figure 1 as a periodic set of parallelograms with side 
lengths T\ and T 2 . The physical periods T[ and T 2 in- 
dicated in figure 1 are defined as the distances between 
crossing points of the line in the nt direction with the 
lines of the periodic structure. These periods can be ex- 



FIG. 3: The light M 2 mass of the scalar field eq.© for the 
Hubble parameter H — 10 _61 tp. The dashed line indicates 
Cr and Ci region for the light mass M2 ~ 2.8 ■ 10 _3 eV and 
the dotted line corresponds to the set of c r and d values for 
Mi « 3 • 10 2 GeV. 



pressed through values of Ti, T 2 and k as follows: 

, . . ImT 2 ReTi - ReT 2 ImTi 

±i = \k\ ; (6) 

ImK ReTi — ReK ImTi 

. , ImTi ReT 2 - ReTi ImT 2 

l l = \k\— . 

ImK KeT 2 — ReK ImT2 

The geometry of periodic structure shows that even 
for comparable values of Ti and T 2 periods the physical 
periods T[ and T 2 of the non-linear oscillator eq. (j3]) can 
differ by ten or more orders of magnitude depending on 
the relative directions of k, Ti, and T 2 . 
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The values T\ and Ti will depend on Hubble param- 
eter H, self-interaction constant A, scale factor a, and 
complex valued constant ci = c r + ia. All values will 
be expressed in Plank units. The value for parameter c\ 
can be estimated from the mass of the fluctuating field 
M<j, s» 3 • 10 _3 eV related to the current value of the cos- 
mological constant. Figures 2 and 3 plot the behavior of 
the light and heavy scalar masses for the current value 
of the Hubble parameter H = lO -61 ^ 1 . It is clear from 
these figures that the heavy mass M\ rs 3 • 10 2 GeV and 
the light mass M 2 ~ 2 — 3 • 10 _3 eV can be found si- 
multaneously for some value of parameter c\. The real 
component c r of this parameter will mostly influence the 
heavy mass and its imaginary part c, will mostly affect 
the light one. We set the values of the constants c r and 
Ci based on the known value of the cosmological vacuum 
energy estimated from H 2 and the wide predicted range 
of supersymmetry breaking ~ 10 2 — 10 3 GeV. Similar- 
plots were also done for H = 10~ 40 tp with analogous 
results. A fundamental and clear physical and mathe- 
matical criteria for the choice of the constant c\ is yet to 
be found. 

The direct calculation of the time dependence of the 
Jacobi sn function for the current value of Hubble pa- 
rameter H = lCT 61 ^ 1 as well as for the H = lO" 40 ^ 1 
for the first milliseconds of the universe exhibit the peri- 
odic behavior of the solution eq. (@| with the periods © ■ 
Figures 4 and 5 demonstrate this dual-periodic behavior 
for the low and high frequency oscillations for the specific 
choice of the parameter c±. The high frequency mode is 
very close to harmonic oscillations (although the second 
derivative clearly shows the deviations from the harmonic 
behavior) . It is extremely hard to plot the low frequency 
processes with a background of a higher amplitude and a 
higher frequency (up to 10-16 orders of magnitude) oscil- 
lations. We used moving average filtration to remove the 
high frequency oscillations. The highly non-harmonic, 
periodical pulse process for H = lO -40 ^ 1 is shown in 
the figure 5. 




FIG. 4: The picture presents the high frequency, heavy 
mass, component of the scalar field Q. The left plot is for 
the Hubble parameter H = 10 — tp 1 . For this case the 
period is T[ « 10 10 'tp and the corresponding mass 
component is Mi = 10" 10 Mp » 10 9 Gev . The right plot is 
for the Hubble value H = 10 _61 tp 1 with the corresponding 
period T[ = 3.2 • 10 1(i tp and corresponding mass 
Mi « 3 • lCT 17 M P « 300 Gev . 
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FIG. 5: The picture presents the low frequency, light mass, 
component of the scalar field Q for H = W~ °tp 1 . For this 

case the period is T 2 ~ 2 ■ 10 22 tp and the corresponding 
mass component is M 2 = 5 • 10~ 21 M P « 50 MeV. The low 
frequency process is highly non-harmonic. 



IV. CONCLUSION 

We were able to find a non-perturbative solution for 
the scalar field with an unusual two-faced behavior. This 
solution has two characteristic mass scales related to the 
Hubble parameter as H 1 / 4 and H 1 ! 2 . Even the possibil- 
ity of a dual- frequency solution with one high mass value 
for a scalar superpartner of N=l SUSY multiplet is itself 
a strong breaking of supersymmetry. It is important to 
notice that these two mass scales are not related to the 
usual energy levels of quantum oscillators and represent 
the fundamental, inner property of the scalar solution 
(HJ). Possibly, in some cases this solution could be re- 
garded as strongly coupled oscillators with two different 
resonant frequencies. Interestingly, the existence of this 
solution with its high and low masses is based on the ap- 
pearance of a tiny term proportional to H 2 in the wave 
equation. This term is related to the third natural scale 
with masses ~ 10 _33 eV connected to the inner frequen- 
cies of the entire universe. The possibility of such scalars 
with the masses quantized in terms of cosmological con- 
stant A = ZH as j^j = n in the framework of extended 
supergravity in de Sitter vacua was discussed in [l3|, [l4[ ■ 

The results presented above leave a few unresolved 
questions. It is important to notice that the deriva- 
tion of the above results was based on a quasi-classical 
rather than quantum approach. Our consideration of a 
static case allows us to replace the non-linear field equa- 
tion with the non-linear oscillatory equation. One of the 
major questions is the observational absence of the light 
scalar particles with the masses 2.5 • 10~ 3 eV Apart from 
the obvious observational difficulties detecting such parti- 
cles there exist additional circumstances that could make 
the detection especially difficult. This light scalar will in- 
teract with its own fluctuations with the same energies 
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as the mass of the particle. In this case an additional 
non-linear dumping will inevitably appear, similar to the 
non-linear attenuation of the waves in a diffuse noisy 
field. Therefore, the life time of this light particle can 
be very short, making it hard to detect. 

Another question is related to the quantum fluctua- 
tions of the high frequency scalar component. Clearly, a 
truly exhaustive answer can be obtained only in a frame- 
work of a supersymmetric quantum field theory in ex- 
panding cosmology. 

The solution described in this paper, where the scalar 
superpartner has such an unusual structure is consider- 
ably novel. If there is a similar structure for the vector 



multiplets is also not clear and can be a subject of future 
research. 

The processes in the early universe during the infla- 
tion are related to fast changes of Hubble parameter as 
well as the scale factor a, and its derivatives a and a 
|9( . The linear term in eq. ([3]) was relatively large at this 
time. Additionally, the sign of this term will depend on 
- . This creates a potential feedback that can induce sta- 
bilizing/ destabilizing effects on the expansion of universe. 

Work of M. Shmakova was supported in part by the 
U.S. Department of Energy under contract number DE- 
AC02-76SF005I5. 
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